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Abstract 

We study the structure of stationary and axisymmetric metrics solving the vacuum Ein- 
stein equations of General Relativity in four and higher dimensions, building on recent 
work in hep-th/0408141. We write the Einstein equations in a new form that naturally 
identifies the sources for such metrics. The sources live in a one-dimensional subspace and 
the entire metric is uniquely determined by them. We study in detail the structure of 
stationary and axisymmetric metrics in four dimensions, and consider as an example the 
sources of the Kerr black hole. 
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1 Introduction 

In recent years there have been a great deal of attention towards research in black holes 

in higher-dimensional space-times. For four-dimensional asymptotically flat space-times 

the Uniqueness Theorems states that only one type of black holes exists for a given set of 

asymptotic charges. In particular, for four-dimensional gravity without matter, the Kerr 

black hole ^ is the unique black hole solution for a given mass and angular momentum 

[H El 12 IS]- Contrary to this, we have learned that for more than four dimensions there are 

generically many available phases of black holes. ^ For Kaluza-Klein spaces, i.e. Minkowski- 

space times a circle, the phase structure is very rich with interesting phase transitions 

between different kinds of black holes, and in some cases even an uncountable number of 

different phases are available For five-dimensional asymptotically flat space-times 

^Notice though the uniqueness theorems of |^ Q |S| for higher-dimensional static black holes in flat 
space-times. 
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without matter we have in addition to the Myers-Perry rotating black hole solution 
also the recently discovered rotating black ring solution of Emparan and Reall \12\? 

To understand the phase structure of black holes in higher dimensions it is important to 
search for new black hole solutions. But the non-linearity of General Relativity makes it in 
general very hard to find new solutions. However, a class of metrics for which the Einstein 
equations simplify considerably are the stationary and axisymmetric metrics. These are 
D-dimensional metrics possessing D — 2 mutually commuting Killing vector fields. Among 
black hole metrics of this type are the Kerr rotating black hole in four dimensions pQ, the 
Myers-Perry rotating black hole in five dimensions [11) and the rotating black ring jl2j . 
In a canonical form for stationary and axisymmetric metrics were found for which the 
Einstein equations effectively reduce to a differential equation on an axisymmetric D — 2 
hy D — 2 matrix field living in three dimensional flat space. ^ The results of generalizes 
those of Papapetrou |15| I16j for stationary and axisymmetric metrics in four dimensions. 
Moreover, the results of also generalizes the work of |17[ ll8j on the Weyl-type metrics, 
which are D-dimensional metrics with D — 2 mutually commuting and orthogonal Killing 
vector fields. Ref. UHl^ is the higher-dimensional generalization of the analysis of Weyl 
jl7j on static and axisymmetric metrics in four dimensions. 

In [Tl] it is furthermore shown that any stationary and axisymmetric solution have a 
certain structure of its sources, called the rod-structure. The rod-structure is connected 
to the three dimensional flat space that the reduced Einstein equations can be seen to 
live in. In this three-dimensional flat space, the sources for the metric consist of thin rods 
lying along a certain axis of the space. To each of these rods can be attached a direction 
in the D — 2 dimensional vector space that the Killing vector fields are spanning. The 
analysis of stationary and axisymmetric solutions in terms of its rod-structure generalizes 
the analysis of UH] for Weyl-type metrics, which again generalizes the analysis of four- 
dimensional Weyl-type metrics (see for example |21j and references therein). 

In this paper we find a new formulation of the reduced form of the Einstein equations 
found in . In the new formulation the sources for a given stationary and axisymmetric 
solution are naturally identified. These sources consist of two D — 2hy D — 2 matrix valued 
functions living on an axis in the above-mentioned three-dimensional flat space that the 
reduced version of the Einstein equations lives in. It is argued that these sources uniquely 
determines the solution. One can therefore say that the problem of finding a stationary 
and axisymmetric solutions is reduced to a one-dimensional problem, i.e. the problem of 
finding the appropriate sources for the solution. 

We examine the general properties of the sources, which are intimately connected to 

the rod-structure of a solution. As an important part of this, we describe how to continue 

^See also the recent review |13| and references therein. 

^Apart from a subclass of metrics that we will not treat in this paper |14|. 

*See Section 4.1 in 19' for a brief review of generalized Weyl solutions. See furthermore 20 for work 
on extending the generalized Weyl solutions of |18| to space-times with a cosmological constant. 
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the sources across rod endpoints. It was conjectured in ^] that a given rod-structure 
corresponds to a unique stationary and axisymmetric metric. The idea being that the 
rod-structure for a solution contains all information about the solution. However, we do 
not find in this paper sufficient constraints to fix the sources in terms of the rod-structure. 
We leave this as a problem for future research. 

For completeness we note here that there have been several recent interesting devel- 
opments in the study of stationary and axisymmetric solutions in higher dimensions. In 
mi 123] a new tool to categorize stationary and axisymmetric solutions by drawing the so- 
called Weyl Card diagrams have been developed. Moreover, solution generating techniques 
for five-dimensional solutions have been explored in |241l25l05] . In particular, in it is 
shown using that Einstein equations are integrable for stationary and axisymmetric 
metrics, and this is subsequently employed in rederiving the Myers- Perry five-dimensional 
rotating black hole. 

In Section [21 we review the results of Jl] that we build on in this paper. In Section 
|31 we find the new formulation for the Einstein equations, and we identify the sources 
for the stationary and axisymmetric metrics. We furthermore examine the properties of 
these sources. In Section IH we prove in the four-dimensional case that stationary and 
axisymmetric solutions are uniquely determined by their sources. We also comment on 
the higher-dimensional cases. In Section we find a way to relate the behavior of the 
sources on each side of a rod endpoint. We prove that we in principle can make a full 
determination of the sources on one side of a rod endpoint by knowing the sources on the 
other side of the endpoint. In SectionElwe consider the asymptotic behavior of the sources 
for asymptotically flat space-times. In SectionjTlwe consider the example of the Kerr black 
hole, flnd its matrix-valued potential and its sources, and we consider the behavior of the 
sources near the endpoints of the rods. In Section |H1 we present our conclusions. 

2 Review of stationary and axisymmetric solutions 

We review in the following some of the most important results of |E] for use in this paper. 
As mentioned in the Introduction the analysis and results of jl4j builds on and generalizes 
results of Refs. ^lElElE]- See the Introduction or ^l] for details on this. 

2.1 Canonical metric and equations of motion 

We consider in this section D-dimensional stationary and axisymmetric solutions of the 
vacuum Einstein equations. These are Z)-dimensional Ricci-flat manifolds with D — 2 com- 
muting linearly independent Killing vector fields V(j), i = 1, ...,D — 2. For such solutions 
we can find coordinates x*, i = l,...,L> — 2, along with r and z, so that 




_d_ 
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for i = 1, .., D — 2, and such that the metric is of the canonical form |14j ^ 

D-2 

ds^ = Gijdx'dx^ + e^^idr"^ + dz"^) , (2.2) 

with 



r = ^|det(G.,)| . (2.3) 
For a given r and z we can view Gij asaD — 2hyD — 2 real symmetric matrix with G*-' 
as its inverse. With this, the vacuum Einstein equations R^i, = for the metric ()2.2() with 
the constraint H2.3|) can be written as 

G-i (^d^ + ^dr + G = {G-^drGf + (G-ia,G)2 , (2.4) 

dri^ = -1 + ^ Tr {{G-^drGf - {G-^d,Gf) , = ^Tr {G-'drG G-'d,G) . (2.5) 

From this one sees that one can find stationary and axisymmetric solutions of the vacuum 
Einstein equations by finding a solution G{r,z) of the equation 1)2. 4() . Given G{r,z) one 
can then subsequently always find a solution z^(r, z) to 1)2. 5|1 since these equations are 
integrable. 

We can make a further formal rewriting of (|2.4|) by recognizing that the derivatives 
respects the symmetries of a flat three-dimensional Euclidean space with metric 

dsl = dr"^ + r'^d-f'^ + dz^ . (2.6) 

Here 7 is an angular coordinate of period 2it.^ Therefore, if we define V to be the gradiant 
in three-dimensional flat Euclidean space, we can write (|2.4|) as 

G-^V^G={G-^VGf . (2.7) 

For use in the paper we elaborate here a bit more on the precise meaning of the formula 
1)2. 7|1 . Using the auxilirary coordinate 7 we can deflne the Cartesian coordinates ui and 
(T2 as 

o"! = rcos7 , CJ2 = rsin7 . (2-8) 

Clearly, we have then that (o"i , cr2 , z) are Cartesian coordinates for the flat three-dimensional 
Euclidean space defined in (|2.6|) since the metric in these coordinates is 

dsl = + d(^2 + dz"^ ■ (2-9) 

We write the Cartesian components of a three-dimensional vector as w = {wi,W2,Wz), 
and the gradiant used in (|2.7|) is then given by 

V = {di,d2,d,) , (2.10) 

where di = d/dai and 82 = d/da2- 

^Two further assumptions are required as premises in the derivation (see [141 for explanation): (i) The 
tensor V^i^j^ Vj2) ' ' ' '^CD'-2)-^'^^(i) vanishes at at least one point of the manifold for a given i = 1, D — 2. 
(ii) det(G'ij) is non-constant on the manifold. 

^It is important to remark that 7 is not an actual physical variable for the solution H2.2|l . but rather an 
auxilirary coordinate that is useful for understanding the structure of Eqs. H2.4|l . 
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2.2 Rod-structure of a solution 



If we consider the condition (|2.3() we see that det(G(0,2;)) = 0. Therefore, the dimension 
of the kernel of G{0, z) must necessarily be greater than or equal to one for any z, i.e. 
dim(ker(G(0,z))) > 1. 

In Jl] we learned that in order to avoid curvature singularities it is a necessary condi- 
tion on a solution that dim(ker(G(0, z))) = 1, except for isolated values of z. We therefore 
restrict ourselves to solutions where this applies. Naming the isolated z- values as ai, a^r, 
we see that the z-axis splits up into the + 1 intervals [— cxD,ai], [ai,a2\-,-.--\aN-i-,ci'N]-, 
[ajv,oo]. For a given stationary and axisymmetric solution we thus have that the z-axis 
at r = is divided into these intervals, called rods. 

Consider now a specific rod [zi, Z2]. From we know that we can find a vector 

" = (2.11) 

such that 

G(0,z)t; = 0, (2.12) 

for i = 1, D — 2 and 2; G [zi, 2:2]- This vector v is called the direction of the rod [zi, 22]- 
We note that if GijV^v^ /r"^ is negative (positive) for r ^ we say the rod [zi,Z2\ 
is time-like (space-like). For a space-like rod [zi,Z2] we clearly have a potential conical 
singularity for r ^ when z £]zi,Z2[- However, if is a coordinate made as a linear 
combination of x', i = 1, — 2, with 

d i d 

drj ^ dx'' ' ^ "^^ 

then we can cure the conical singularity at the rod by requiring the coordinate ry to have 
the period 



A?? = 27rlim^^/^— - , (2.14) 

for z e [zi,Z2]. 



3 New formulation and identification of sources 

In this section we rewrite the vacuum Einstein equations in a way that gives a natural 
identification of the sources for a given stationary and axisymmetric solution. The purpose 
of identifying such sources is that one can hope to reduce the problem of finding new 
solutions to the problem of finding the sources for new solutions. We comment on the 
general philosophy behind the introduction of these sources in the end of this section. 
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3.1 New formulation 

We begin by defining the field C(r, z) by 

C = G-^VG. (3.1) 

Tlius, Cr = G-^drG and = G-^d,G. With this definition, the equation ^Jii for 
G{r, z) is equivalent to the equation V • C = for r > 0. However, for r = we can have 
sources for C. Therefore, we require C{r,z) to obey the equation 

V -C = Att6'^{(t)p{z) , (3.2) 

i.e. that r~^dr{rCr) + dzCz = 'i'Kd'^{a)p{z). Here d'^{a) = d{ai)6{a2), where ai and ct2 are 
the Cartesian coordinates defined in (|2.8|) . Thus, the delta- function 6'^ (a) expresses that 
we have sources for C{r, z) at r = 0. The D — 2 hy D — 2 matrix- valued function p{z) in 
H3.2|) then parameterizes the sources for C(r, z) at r = 0. 

Interestingly, the equation H3.2|) is now a linear first order differential equation which is 
equivalent to the non-linear second order equation l|2.7|l . However, the price of introducing 
C is that C in addition should obey the non-linear equation^ 

V xC + C xC = , (3.3) 

i.e. that drC^ + d^Cr + [C^, C^] = 0. The equation is found from V x (G^^VG) = 
VG-i X VG = -G-^VGG-^ x VG = -(G-^VG) x (G-^VG). Note that there are no 
source terms to the equation (|3.3jl so it also holds for r = 0. In conclusion, we have 
exchanged the non-linear equation second order equation (|2.7|) with the two first order 
equations ()3.2() and (|3.3() . 

We see from Eq. ()3.2|) that Cij{r,z) resembles an electric field with Pij{z) being the 
charge density at r = 0. However, the C x C term of Eq. (|3.3j) introduces terms that 
mix the different components in such a way that C{r,z) cannot be obtained with linear 
superposition of fields. 

The constraint 1)2. 3|) on the determinant of G is equivalent to demanding 

Tr(G) = 2Vlogr , (3.4) 

i.e. that Tr(Gr) = 2/r and Tr(G^) = 0. Therefore, if we take the trace of (|3.2)) we get, 
using (|3.4() and V^logr = 2-K5'^{a), the important identity 

Tr(/>) = 1 , (3.5) 

which holds for all z. Eq. 1)3. 5(1 is a generalization of the requirement on the Generalized 

Weyl solutions of jl8j that the total rod density is constant for all z. As shown in ^1], we 

'^Note our conventions for the cross product is such that the epsilon symbol has e^^^ = 1 which in 
{r,"f,z) coordinates means that e''^^ = 1. 
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get from this and demanding absence of singularities for r ^ that there is precisely one 
rod present for any given value of z (except for isolated values of z). 

Note that we can express the equations (|2.5() for i/(r, z) in terms of C(r, z) as 

dr^ = -7^ + ^ 1V(C2 - d) , 8,1^=^- Tr(aC,) . (3.6) 
2r 8 4 

These equations ensure that we can find z^(r, z) from C(r, z) alone. 

It is interesting to consider how many different G{r, z) that corresponds to the same 

C(r, z). The answer is provided by the following lemma: 

Lemma 3.1 Let Gi{r,z) and G2{r,z) be given such that G^^VGi = G2^VG2, with Gi 
and G2 being invertible for r ^ 0. Then we have that G2Gi^ is a constant matrix. 
Proof: Define M = ^2^5"^ Then 

G-^VGa = (MGiy^ViMGi) = G^^M-i((VM)Gi + MVGi) 
= G^iM-i(VM)Gi + Gf^VGi . 

We see from this that VM = for r 7^ 0. This means that M is a constant matrix. □ 

The above lemma means that given C(r, z) we can determine the corresponding G{r, z) 
up to a constant matrix, hence if we for example require certain asymptotic boundary 
conditions on G(r, z) then we can determine G(r, z) uniquely for a given C(r, z). 

3.2 Defining a potential for C 

We define the D — 2 hy D — 2 matrix- valued function A(r,z) by 

Cr = --d,A , C, = -drA . (3.7) 
r r 

This definition is meaningful for r > since r~^dr{rGr)+dzCz = implies the integrability 
condition drdzA = dzdj-A. Hence, for a given C we can find A such that 1)3. 7() is fulfilled 
for r > 0. The value of A for r = is then found by demanding continuity of ^4 at r = 0. 

Clearly, A(r,z) is a potential for G{r,z). If we define the field A(r,z) by Aj.{r,z) = 
Az{r, z) = and A^{r, z) = A{r, z) we can furthermore write that C = V x A. Therefore, 
the potential A{r,z) can be seen as a component of the vector potential A{r,z).^ Note 
also that the definition of A from C is ambiguous since if one uses the vector potential 
A' = A + ^a for a given a(r, 7, z) one gets the same C. However, demanding A'^ = A'^ = 
and that d-yA' = we get that the definition 1)3 .71) for A{r, z) is only ambiguous up to an 
additive constant. 

We can now use Gauss' law for C to derive an important identity. Consider the 
cylindrical volume V = {r < r^, zi < z < Z2}. Then we have 

— / V • CdV = — [ n- CdS 

27r Jy 27T Jqv , , 

rro f-ro i-Z2 

= - / drA\z=z^dr + / drA\z=z2dr - / dzA\r=r^dz = ^(0, Zi) - A(0, Z2) . 

J r=0 J r=0 J z=z\ 



Note that in Cartesian coordinates A = (j4i,^2, Az) — (j49i7, j4927, 0). 
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Using (|3.2|1 we then get^ 

p{z) = -]^d,A{Q,z) . (3.9) 



Therefore, we get from p.7p that 



p{z) = \\imrCr . (3.10) 



We see then that 2p = hm^— >o ?'G^^(?r G which means that one can compute p from the 
metric directly. 

We have thus solved ()3.2() by introducing the potential A{r, z) via ()3.7() and demanding 
(|3.9|) . However, we still have (|3.3j) which for A{r,z) becomes 

rdr {^d' A^ + dlA + ^[drA, d^A] = . (3.11) 

Thus, ()2.7() for G{r,z) have now been translated into (|3.11() for A{r,z). The constraint 
(|2.3)) . which is equivalent to (|3.4)) . translates for the potential A{r,z) to 

Tr(^) = -2z , (3.12) 

up to an additive constant. 



3.3 Behavior of a solution near r = 

Equation (|3.1U|) shows that p{z) characterizes the behavior of Cr{r, z) for r — > 0. However, 
it will be clear in the following that we in addition need to take into account the behavior 
of Cz{r, z) for r — > 0. We define therefore the D — 2 by -D — 2 matrix- valued function A.{z) 

by 

A(z) = C,(0,z) . (3.13) 
We comment on the finiteness of A{z) below. From (|3.4|) it is easily seen that 

Tr(A) = . (3.14) 

Moreover, from considering the r — > limit of H3.3|) . we get 

p' = [p,A]. (3.15) 

Therefore, if one knows A{z) it is possible to find p{z), at least up to an additive constant 
matrix. 

Consider now a rod [zi,Z2], as defined in Section 12.21 Without loss of generality we 

consider the rod [zi,Z2] to have direction v = (1,0, ...,0). We consider what happens for 

general v below. Using the analysis of the metric near a rod in we see then that for 

''The identity (gSJ can also be derived using the Cartesian coordinates l|2.8|l . Here one uses [di,d2]'y ~ 
27r5^((T). Note that just as for the Dirac monopole in electrodynamics one encounters the Dirac string at 
r — 0. The integral version of the derivation 13.811 instead avoids this issue. 
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r ^ and z ^]zi,Z2[ we have that Gu = 0{r'^), Gu = Gn = 0{r^) and Gij = 0(1) with 
i,j = 2,...,D - 2, while for the inverse we have that = 0{r~^), G^' = G'^ = 0(1) 
and G*-' = 0(1) with i,j = 2,...,D — 2. Moreover, writing Gu ~ ite'^^^V^, we have 
that G^^ ~ ^e^H^)r-'^ for r ^ and z €z]zi,Z2[- From l\3.1\\ this gives then for r — > 
and z £]zi,Z2[ that Cr,ii — l/r, Gr,ii = 0{r^^), Cr,ii = 0(r) and Gr,ij = 0(r) with 
i,j = 2,...,D- 2, while C,,ii ~ A'(z), = 0(1), = 0{r^) and C,,^^ = 0(1) with 
i,j = 2,...,D-2. 

Using then (|3.10() and (|3.13() we conclude that both p{z) and A(z) are well-defined 
and smooth for z £]zi,Z2[. Moreover, for p{z) we get that pu = 1 and that pij = with 
i = 2, D — 2 and j = 1, ...,D — 2 for z €]zi, Z2[, while for A{z) we get that An = A'(z) 
and Aji = with i = 2, ...,D — 2 for z €]zi,Z2[- Note that we can infer from this that 
Tr(pA) = A'. This is used below. 

Consider now instead a given rod [zi, Z2] of direction v, with v being arbitrary. Clearly, 
we can find X{z) such that v^Gv = ■ii\{z)r'^v^ v. Using now the above analysis along with 
Appendix^ we can conclude for p{z) with z £]zi,Z2[ that 

pv = v , Tt{p) = 1 , p^ = p , (3.16) 

and that there exists linearly independent vectors wi, wd-s independent of z such that 

p'^Wi = 0,1 = 1,...,D- 3 . (3.17) 

We see from 1)3. 16() that v is an eigenvector for p with eigenvalue one. Concerning the 
properties involving A, we have for z ^]zi, Z2[ that 

Av = X'v , Tr(/>A) = A' . (3.18) 

Thus, V is also an eigenvector for A(z). If we consider i/(r, z) for r ^ 0, we see from ()2.5() 
that dri' — —(1 — Tr(/9^))/(2r). Using now (|3.16jl we see that Tr(/9^) = 1 and hence that 
z^(0, z) is well-defined for z £]zi,Z2[- If we instead consider d^iy for r — > we get from 
and (jTTK|) that 

a,KO,^) = ^TV(pA) = ^A'(z) , (3.19) 

for z G]zi,Z2[- Integrating this, we get 2z^(0, z) = X{z) + constant. This ensures that 
for space-like rods we can remove the potential conical singularity by choosing a suitable 
periodicity of the corresponding coordinate, while for time-like rods it ensures the temper- 
ature of the corresponding horizon is constant on the horizon. We see furthermore from 
(|3.19j) that we can compute z^(0, z) for all z by knowing p{z) and A(z), and demanding 
continuity of z^(0, z). 

General philosophy 

The general philosophy of introducing the sources p{z) and A(z) is that all information 
about the solution is stored in these two matrix-valued function. This is shown explicitly 
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in the four-dimensional case in Section^ Thus, we have reduced the relevant physics to 
these two matrix-valued functions living on a one-dimensional subspace. In particular, all 
the physical quantities like the mass, angular momenta, temperature, entropy, quadrupole 
moment, etc. can be obtained from p{z) and A(z) directly. We show this explicitly for 
the most important asymptotic quantities in Sectional 

The idea is moreover that the sources p{z) and A(z) should be uniquely determined 
by the rod-structure, i.e. by specifying the division of the z-axis into separate intervals, 
along with specifying the direction of these rods in the linear vector space of Killing 
vector fields. This would be a generalization of the solution generating technique for Weyl 
solutions j^E]. In this paper we do not go all the way to determine p{z) and A{z) 
from the rod-structure, but in addition to the steps taken in this section we take a further 
important step in this direction in Section [SI 

4 On uniqueness of A{r, z) given p{z) and k.{z) 

The purpose of this section is to argue that stationary and axisymmetric solutions of the 
vacuum Einstein equations are uniquely determined by the sources p{z) and h.{z) intro- 
duced in Section 131 We address this here by arguing that the potential A{r, z) introduced 
in Section 131 can be uniquely determined by p{z) and A(z). 

We prove in this section in detail that A{r, z) is uniquely determined by the two matrix- 
valued functions p{z) and A(2:) in four dimensions. This means that we can reduce the 
question of finding a solution for A{r, z) to finding the two functions p{z) and A(z) that 
lives on a one-dimensional subspace. We comment on the proof for arbitrary dimensions 
in the end. 

The general idea is to expand A{r, z) around r = 0. The expansion is 

oo 

A(r,z)=^>l„(z)r2". (4.1) 

n=0 

Since the behavior of A{r, z) is singular around rod endpoints we assume zi < z < Z2 
where [zi, Z2] is one of the rods of the solution. We get then from (|3.11|1 

^^; = 2K,Ai], A2 + \[A2,A'^] = -\A'i-\[Ar,A'^], (4.2) 

^n+l + ^[^n+r, A^] = - • (4.3) 

Note that A'^ = —2p and Ai = ^A. Therefore, the basic idea in the following is to show 
that we can find all An for n > 2 from A'q and Ai. 

We now restrict ourselves to four dimensions. We can assume without loss of generality 
that we have a rod [zi, Z2] in the direction 
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> /3=(::i , 7=(::i . (4.5) 




Define 

^0 0^ 

These three matrices have the properties 

= a , = , 7^ = , a/3 = P , /3a = , a^y = , -ya = 7 , 

(4.6) 

Pl = a, -ff3 = I-a, [a,/3]=/3, [0,7] = -7 , [/?,7]=2a-/. 

From the properties of p of Section r3.3l we see that it should have the form p(z) = a + h{z)l3 
where h{z) is a function. Moreover, the most general form of A is A(z) = \'{z){2a — I) + 
{h'{z) + 2h{z)\'{z))f3 where \{z) is defined in Section lOl Note that ^0 = -2a - 2/i/3. 
Since Tr(^„) = for n > 1 we can write 

An = an (2a - /) + 6„/3 + c„7 , (4.7) 

for n > 1 and zi < z < Z2- Note that ai = ^q', bi = ^h' + hq' and ci = 0. From 1)4.31) we 
get then for n > 1 



4n(n + l)a„+i + 4(n + l)/iCn+i + + 2^(n - A: + l)(6„_fc+i4 - Cn-k+ib'k) = , 

fc=i 

n 

4(n + l)^6„+i - 8(n + l)/ia„+i + 6^ + 4^(n - /c + l)(a„_fc+i6'fc - 6„_fc+i4) = , 

k=l 

n 

4(n - l)(n + l)c„+i + c'n + 4^(n - A; + 1)(c„_a..+i4 " On-fc+iCfc) = . 

k=l 

(4.8) 

For n > 2 we see that these three equations uniquely determine a„+i, bn+i and c„+i from 
Ofc, 6fc and for k < n. However, for n = 1 these equations do not constrain C2{z). 
Therefore C2{z) can be an arbitrary function, according to these equations. 

In order to fix C2{z) we should impose that we want A{r, z) to be a potential for C(r, z) 
such that C = G^^VG with G^ = G. Then we know from Appendix ^ that h{z) and 
\'{z) can be chosen freely. On the other hand, from the expansion of G{r,z) in Appendix 
IbI h{z) and \'{z) are also the only free functions. Therefore, C2{z) must be determined in 
terms of h and A'. Indeed from (|B.12() we get 

C2{z) = --/'h' . (4.9) 

So, given h{z) and X'{z), we must choose C2{z) as ()4.9)) in order for A(r, z) to be a potential 
for which a symmetric G(r, z) exists. Said in another way, imposing 1)3.11(1 on A{r, z) 
ensures the existence of a G(r, z) such that C = G^^VG, but imposing in addition that 
G^ = G means that we also have to fix C2{z). 
We can formulate the above theorem: 

Theorem 4.1 Let a rod \z\^Z2\ be given. Given p{z) and A{z) for zi < z < Z2 there exists 
precisely one A{r,z) that solves ()3.11|) and that corresponds to a symmetric G{r,z), for 

Zl < Z < Z2. 
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To show that this theorem works in arbitrary dimensions one should generahze the 
procedure above. One can of course without loss of generalitity choose a rod of direction 
(1,0, ...,0) generalizing the four dimensional case. The proof basically rest on having 
that An^i for a given n is determined by Eq. (|4.3|) . If we take into account the general 
properties of p from Section ESI we see that this is only problematic for n = 1. For n = 1 
we have that only the top entry of the first column of the matrix A2 — [A2, p] is non-zero, 
no matter what A2 is. Therefore, given p and A we can choose D — 3 additional arbitrary 
functions and still obey ()3.11() . generalizing the arbitrariness of C2{z) for four dimensions. 
We expect these D — 3 free functions to be determined by the = G condition, just as 
in four dimensions. 

5 Smoothness conditions at rod endpoints 

In this section we take a further step in understanding how to obtain the sources p{z) and 
A(z) from specifying the rod-structure for a solution. What we address in this section is 
how to relate the sources on both sides of a rod endpoint. Using this one can find the 
sources on one side of a rod endpoint from knowing the sources on the other side of the 
rod endpoint. 

We can address the question of how to connect the sources across rod endpoints in a 
mathematical way as follows. We have seen that p{z) and A(z) are smooth matrix valued 
functions away from endpoints of rods. But, this is not true at the rod endpoints. This 
is directly related to the fact that the potential A{r, z) is not smooth near rod endpoints. 
However, we want the solution as a whole to be regular, thus we need a criteria for when 
a solution is well-behaved near a rod endpoint. This is the subject of this section. 

Our basic requirement for smoothness of a solution is: 

• Let be given. Define the coordinates 



We then require of a solution that A{p, q) should be smooth in a neighborhood of 



Away from the endpoints, this reduce to usual requirement of smoothness of A{r, z) near 
r = 0, which already is contained in the previous sections. However, if z-t is an endpoint, 
i.e. if we suppose that we are considering a solution with two different rods [zi,^*] and 
[z*, Z2] so that z = divides the two rods, then the above requirement becomes non-trivial. 

In the following we use the above smoothness requirement to analyze the behavior of 
p{z) and A{z) across endpoints. Let now be given, and consider the (p, q) coordinates 
as defined in (|5.1|) . We define 



p = z — z^ 





{p,q) = (0,0). 



rp{m,n) ^ Qrngn^ 



'p '^<?^l(p,g)=(0,0) 5 



(5.2) 
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as notation for the derivatives of 74(p, q), with T^'^'^) =74(0,0). Using that 2p = —dpA\r=o— 
sgn{p)dqA\r=o and A = "^95^1^=0, we find that the expansions of p and A around z = 
are given by 



^ = ^ E It E (- sgn(p)'r(-^-^.o _ sgn(p)^+^r(-M-^)) , (5.3) 

M = E It E Ju^ sgn(,)'+^T(-M+^) . (5.4) 
k=o ' 1=0 '' 

If we consider T'^'^'^^ we see from (|5.3j) and (|5.4|) that 

hm p{z) — lim p{z) = hm {z — z^,)K{z) = — hm {z — z*)A(z) . (5-5) 

2— >2:7 Z^zt Z^z't Z^ZZ 

Considering instead T(1'0) and T^^.o) _ ^(0,2)^ 

Mm [2p(2) + (z - z*)A(z)] = lim [2/9(z) + (z - z,)A(z)] , (5.6) 
lim [p'(z) + ((z - z,)A{z)y] = lim [p'(z) + {{z - z*)A(z))'] . (5.7) 

Eqs. ()5.5() - H5.7p relates the behavior of p{z) and A{z) for z ^ to their behavior for 
z ^ z~ . This means that knowing p{z) and A{z) for z > z^, we can get information on 
p{z) and A(z) for z < z^,. The question is now whether one can fully determine p(z) and 
A(z) for z < z*. This is in fact possible though one needs to use the extra constraints on 
coming from expanding 

d^pA + d^A + ^dpd.A = ^ [dpA, dgA] , (5.8) 

around {p,q) = (0,0). Eq. 1)5. 8p is Eq. ()3.11|) in {p,q) coordinates. The first few terms of 
this expansion give 

2^(2,0) _^ 2^(0,2) ^ [2^(1,0)^ 2^(1,2)] _^ [2^(1,1)^ 2^(0,1)] . ^^'^'^ 

That we can connect p{z) and A(z) across a rod endpoint is formulated in the following 
theorem: 

Theorem 5.1 Let two rods [zi,z^,] and [z^,,Z2] be given. Assume that we know p{z) and 
A(z) for zi < z < z^,. Then we can determine p{z) and A(z) for z* < z < Z2 uniquely by 
the smoothness condition on A{p,q) formulated above with {p,q) given in H5.1I1 . 

We can prove this as follows. Since we know p{z) and A(z) for zi < z < z* we can determine 
A(r,z) for zi < z < z* using Theorem 14.11 ^*^ Employing the coordinate transformation 

We assume here the validity of Theorem 14.11 for all dimensions, although we strictly speaking only 
have proven Theorem 14. II in detail in four dimensions. 
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(|5.1|) we can then find A{p, q) for p <Q and g in a certain range. Since A{p^ q) is assumed 
to be smooth in {p, q) = (0, 0) we see that all the derivatives of A{p, q) at {p, q) = (0, 0) 
are the same whether one approaches from p < or p > 0. Thus, we can determine T^"^'") 
as defined by ()5.2|) from our knowledge of p{z) and A(z) for zi < z < z^. Using now ()5.3|) 
and (|5.4j) we can find the values of the derivatives of p{z) and {z — z^,)h.{z) as z ^ z^ to 
all orders. This then uniquely determines p{z) and A(z) for z^ < z < z^- 



6 Asymptotically flat space-times 

In this section we consider solutions that asymptote to four- and five-dimensional Minkowski 
spaces. We work in units where Newtons constant is set to one. 

An important purpose of this section is to show that one can find the asymptotic 
quantities like mass and angular momenta directly from the sources and A(z). This 
is an alternative way of viewing the statement that all the information about a solution 
is contained in the sources /3(z) and A(2;). 



Four-dimensional asymptotically flat space-times 

From ni we know the asymptotic behavior of four dimensional solutions asymptoting to 
Minkowski space. Using this, we get the following asymptotic behavior for A(r, z) 



A 



2Mz . 2Jz{3r'^ + 2z^) . 

11 — / ^ , ) Al2 — , g oxo/o , A2I 



2Jz 



(r2 



i)3/2 ' 



(6.1) 



in the asymptotic region \/r^~+^ 00 with z/^/r^~+~z^ finite. Here M is the mass and 
J is the angular momentum. Combining this with Eq. (|3.8I) and Eq. 1)3.2(1 we get 



(6.2) 



dzpn{z) = 2M , / dzpi2{z) = -4 J , / dzp2i{z) = 

1 J —00 J —00 

Using Eq. ()3.1U|) we get that p{z) and A(^) asymptotically behave as 

/ 0\ \ 



2J 



A(z) = sgn(z) 



/ 2M 
6J 2M 



(6.3) 



for z zbco. Note that the entries with zeroes are zero to all orders in 1/z as can be seen 
from the properties of p and A found in Section 13.31 
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Five-dimensional asymptotically flat space-times 



We know from [Ilj the asymptotics of solutions that asymptote to five-dimensional Minkowski 
space. Using these results we find the following asymptotic behavior of A{r, z) 



AM 



111 



Svr _j_ -,2 



A 



J I \/r2 + z2 + 



J2 \/r2 + 



21 



TT 



j,2 _|_ ^2 



i31 



^ 2 Ji 2zVr2 + z2 + r2 

^12 = — o ) ^22 



TT 



y.2 _|_ ^2 



+ + z) + 



TT 

2(M + r/) 



J.2 _j_ ^2 



2 J2 2zVr2 + z2 _ r2 



13 



TT 



j,2 _|_ ^2 



A 



33 



+ z^ 



z + 



Svr 

2(M - r/) z 



\/r2 + z^ 



23 



r2 + 2z2 + 2z\/r2 + z2 
C T^. ^KTl^ ' ^32 — -C 



(r2 + z2)3/2 



Svr + z2 

r2 + 2z2 - 2zVr2 + z^ 

(^^2 _|_ 2-2^3/2 ' 



(6.4) 



in the asymptotic region \/r2 + z^ ^ 00 with z/\/r2~+^ finite. Here M is the mass 
and Ji and J2 are the angular momenta. See for comments on r/ and Using now 
Eq. together with Eq. and Eq. (|^ we get 



/ dzpii{z) = ——, I dzpu{z) = , / dzpis{z) = 

J-00 -371" J-oo TT 



(iz/92l(z) 



dzp3i{z) 



dzp2z{z) 



dzp32iz) 



. 



(6.5) 



Finally, we see from Eq. 1)6. 4p using Eq. (|3.1UI) that to leading order 





( 








\ 




/o 


\ 




piz) = 




_Ji_ 


1 


2C 


, A(z) = 





1 

— 


+ 




vrz^ 


Z2 




z 

1 

- 

z ' 














J 









for z ^ cxD and 



4M 







37rz2 
3Ji 2(M + r?) 



vrz'" 

A 

vrz^ 



37rz2 




8J< 



vrz" 
8C 

^3 



2(M - 77) 
37rz2 / 



(6.6) 



for z 




J2 2C 



7rz2 z2 



1 



\ 


/O \ 




, A(z) = 


1 

0-0 

z 

(00-lJ 


+ 


/ 





4M 


8J1 


37rz2 


7rz2 


A 


2(M + r?) 


vrz^ 


37rz2 


3J2 


8C 


TTZ^ 


Z3 








37rz2 / 



(6.7) 



-00. 



7 Example: The Kerr black hole 

To illustrate the methods of this paper to analyze stationary and axisymmetric solutions, 
we consider in detail the Kerr black hole pp. 
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The Kerr metric is formulated most easily in the prolate spherical coordinates (x, y) 
defined as 



r = a\J {x^ — 1)(1 — y2) , z = axy , x > 1 , —l<y<l, (7.1 
where a > is a number. The metric is then given as 

-X -XA \ 



G 



-XA X-Hr^ 



X 



cos^ A + sin^ A — 1 



A 



- A^X^) ) ' 

2a tan A(l - y^) (1 + x cos A) 



(7.2) 
(7.3) 



(1 + X cos A)^ + y2 sin^ A ' x^ cos^ A + sin^ A — 1 

We now find the potential A for the Kerr metric. For this we record that in prolate 
spherical coordinates C = V x ^ becomes 



75 —dyA , Cy — — 7T-dxA . 



a{x^-l) ' ' ' a(l-y2) 
Using this we get the following potential for the Kerr black hole 



(7.4) 



A 



11 



2y sin A 



2ay((l + X cos A)^ + sin^ A) 
cos A((l + X cos A)2 + y2 gin^ A) ' (1 + X cos A)2 + y2 sin^ A 

A12 = Ga^ y I 1 



y 



sin A ^ 2a^ 1/(1 — y^)^ sin Atan^ A 



3 / cos^ A (1 + X cos A)2 + y2 sin^ A 

From we have that the rod-structure of the Kerr black hole consist of a space-like 
rod [—00, —a] in the direction (0,1), a time-like rod [—«,«] in the direction (1,0) and a 
space- like rod [q,oo] in the direction (0, 1). Here is given as 

sin A cos A 



n 



2a(l + cosA) 



(7.6) 



We now find the source distributions p{z) and A(z) for the Kerr metric. For z < —a we 
have 



A 



/ 

-(.1, 

-A' \ 
h' + 2h\' A' / 



Hz) 



X'{z) 



2 sin A cos A (l — ^ cos A) 



a 



1 - f cos a) ^ + sin^ A 



2 ' 



2 




1 cos A) ^ - sin^ A 




a cos A ^1 — ^ 




(l-f cosA)^ + sin2A 



(7.7) 



This corresponds to the [—00, —a] rod. For z > a we have 

2 sin A cos A (l + ^ cos A) 




h 1 



Hz) 



a 



1 + f cos A) ^ + sin^ A 



2 ' 



A 



h' + 2hX' X 



2 


;(!+ 


z_ 

a 


COS a) ^ 


— sin^ A 




a cos A ^1 




) 


(1 + ^ 


cos a) ^ + sin^ A 



(7.8) 
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This corresponds to the [a, oo] rod. Finahy for —a < z < a we have 

\-^h h\ \ sinA((l + cosA)2-||sin- 



S72a (l + cosA)2 + 4sin^A 



A 



-J^/i' + (1 - 2Jl/i)A' h' + 2h\' \ (7,9) 

+ 2Q(1 - f]/i)A' rj/i' - (1 - 20/i)A'y ' 
4(1 + cos A)z 



a2(i_£|) (l + cosA)2 + 5sm2A 



A'(z) 



corresponding to the [—a, a] rod. 

One can now easily check that p{z) and A(z) given above for the Kerr black hole 
obeys the properties found in Section l3^ It is also interesting to consider the smoothness 
conditions of Section [21 on the endpoints of the rods. Consider the endpoint z = a. Using 
(|7.9|) we compute 



hm p=\ \ , hm z — a A = 

r r ' n N a ^/^ 1 - 2 COS A / 1 \ 

^hm(p+((.-a)A)) = ^^(^^^ -ij 

Using ()7.8p we compute 

( ^ ^\ V , NA 1 

iim p = \ , iim (z — a)A = 

II I, N AN/N 1 - 2cosA / 1 \ 

hm (p' + ((z - a)A)') = J 

2^a+ 2a \ 217 — 1 / 



(7.10) 



(7.11) 



From (|7.1Up and (|7.1ip one can now easily check explicitly that 1)5. 5^ - 1)5. 7(1 are obeyed. 
Similarly one can check the smoothness conditions involving higher derivatives of p and 
A. 



8 Conclusions 

In this paper we have made a new formulation of the Einstein equations for stationary and 
axisymmetric metrics. This was done in Section |21 by introducing the field C = G^^VG 
and its vector potential A given by C = V x A. This enabled us to naturally identify 
the sources p{z) and A(z) for a given solution. As argued in Section |1] a solution should 
be completely determined once p{z) and A(z) are given. Hence, we have reduced the 
problem of finding stationary and axisymmetric solution to the problem of specifying the 
sources p{z) and A(z). This means that we have reduced the problem of finding solutions 
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to a one-dimensional problem, since p{z) and A(z) are matrix-valued functions of just one 
variable. 

That all information about a solution is contained in the sources p{z) and A.{z) is also 
confirmed by the fact that one can extract all physical quantities from these. In particu- 
lar, in Section El we consider how to find the most important asymptotically measurable 
quantities for four- and five-dimensional asymptotically flat solutions. 

The remaining problem is to find an effective method to determine p{z) and A{z) 
solely from information about the rod-structure of the solution. A step towards this is 
provided in Section |31 in which we analyze the constraints on p{z) and A(z) coming from 
the r — > limit of the Einstein equations, and from imposing a particular rod-structure. 
These constraints significantly reduce the freedom of choice for p{z) and A{z), in particular 
for the four-dimensional case the left-over freedom is in terms of a single function of z. A 
further step towards restricting p{z) and A(z) from the imposed rod-structure is taken in 
Section 13 where we analyze the conditions on p{z) and A{z) when crossing a rod endpoint. 
However, we do not find a sufficient amount of restrictions on p(z) and A(2:) to determine 
them completely from the given rod-structure. We leave therefore this as an open problem 
for future research. 

In Section [7| we applied our considerations to the example of the Kerr black hole. We 
found the potential A{r, z) for the Kerr black hole and from that the sources p{z) and 
A{z). Moreover, we checked successfully that the conditions of Section |S1 on the specific 
sources when crossing rod endpoints are obeyed. 

In conclusion, we have accomplished several steps towards the goal of this paper, 
which is to find a general method of finding stationary and axisymmetric metrics given 
a specific rod-structure. However, we are still lacking a complete determination of the 
sources. Despite this, it is clear from the uniqueness theorems in four dimensions that 
one for example should be able to determine completely the sources of the Kerr black 
hole. On the other hand, we do not now of any direct construction of the Kerr black 
hole. Indeed, if one determines p{z) and A{z) directly from the rod-structure of the Kerr 
black hole it will be the first example of such a direct construction. One can therefore 
say that by devising a method to find the sources solely from the rod-structure one would 
be obtaining a generalization of the uniqueness theorems in four dimensions to stationary 
and axisymmetric solutions in higher dimensions. 

Clearly, the problem of finding the sources from the rod-structure is interesting also 
with respect to the search for new black hole solutions, in particular the five-dimensional 
asymptotically flat black ring with two angular momenta, generalizing ^^1) and non-zero 
angular momenta solutions with Kaluza-Klein bubbles and event horizons TO^. 
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A Rigid transformations 

Consider a coordinate transformation 

x' = ^MijX^ . (A.l) 
i 

Here M can be any (constant) invertible matrix. We see then that 

G = M-^GM~'^ , C = MCM-^ , A = MAM'^ . (A.2) 
Moreover, we have 

3 3 

With a slight abuse of notation we write the first equation of ()A.3|) as 

V = Mv , (A.4) 

i.e. in this equation we read v and v as vectors in instead of vector fields in 

the (-D — 2)-dimensional linear space spanned by the Killing vectors (clearly v^d/dx^ is 
not affected by the transformation (|A.H) ). Thus, the point of the above is that by a 
transformation ()A.1|) we can transform a given non-zero vector v into any non-zero vector 
that we want. We just need to take into account ()A.2|1 . 



B Analysis of EOMs in four dimensions 

In this appendix we consider the expansion of the EOMs (|2.4() for G around r = in four 
dimensions. In doing this, we show that a solution of the EOMs is completely determined 
by two functions. We use this in Section^ 
We parameterize G as 

^=[-XAX-u/-A'X')] • ^^-'^ 



The EOMs (|^ in terms of X and A are 



X{d^ + ^dr + d^^ X = (drX)^ + {d^Xf - ^ ({drAf + {O^Af) , (B.2) 

a, (^drA\ + d, i^d,A\ = . (B.3) 
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(B.6) 



We consider now the behavior of a solution for r small (i.e. near r = 0) and for z away 
from the endpoints of any rod. Then we can expand X and A as 

oo oo 

X(r, z)=Y, /„(z)r2" , i(r, z) = ^ g^^zV^ . (B.4) 

n=0 n=0 

Considering then the leading part of (|B.2|) we get that either /o(-z) is zero, or g^iz) is 
constant. 

We consider first the case in which (70(2) is constant. We write 

/o(z) = e'^(^) , 5o(^) = -a, = -M^)e2^(^) . (B.5) 

The first few terms in X and A then takes the form 

X = e-^(^) + r'h^z) + rV2(^) + O(r') , 
i = -a - r2/i(z)e2^(^) + r^g^^z) + r653(z) + ©(r^) . 
One can now find /i(z) and 52 (^) from the EOMs (|B.2|) - ()B.3|1 : 

/i = ^e-^A" - e^}? , <72 = le^\h" + 4/iA" + 2/i'A') - e^^h^ . (B.7) 
4 8 

By going order by order one can similar find /2(-z) and 53(2), and then fsiz) and 54 (-z), 
and so on. Therefore, once h{z) and X{z) are given, X(r, z) and A{r, z) are uniquely 
determined. We can now find p and A: 

_ f ah a{l - ah)\ ^_ ( a{h' + 2h\') - A' -a'^{h' + 2h\') + 2aA' \ 

1-a/i y' ~y h' + 2hX' -a{h' + 2/iA') + A' J ' ^ ' ^ 

One sees from this that if one knows both p and A then X{r, z) and A[r, z) are uniquely 
determined. (Notice though the constant part of \{z). However, that is not essential since 
that can be scaled to whatever one wants). 

We consider now the case for which fQ{z) = 0. We write 

h{z) = , h{z)= e^(^) , g^iz) = h{z) . (B.9) 



We then get 
Writing 

we also compute that 




, . ,„ , . fx' h' + 2hX'\ 



C,(r, z) = K{z) + Ar^A2(z) + O(r^) . (B.ll) 

^2,21 = --/"-h' . (B.12) 
This is used in Section |11 One can solve for f2{z) and gi{z) from the EOMs (|B.2() - (|B.3() : 

h{z) = -ie^A" , 91 (z) = -h" - \h'X' . (B.13) 

By going order by order one can similar find f^{z) and 52 (-2), and then fA{z) and 93 (z), 
and so on. Therefore, once h{z) and X{z) are given, X(r, z) and A{r, z) are uniquely 
determined. 
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